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ABSTRACT 

^ ! We consider the amplification and transport of a magnetic field in the 

collapsed core of a massive star, including both the region between the 
neutrinosphere and the shock, and the central, opaque core. An analytical 
argument explains why rapid convective overturns persist within a newly formed 
■^j- ! neutron star for roughly 10 seconds (> 10 3 overturns), consistent with recent 



numerical models. A dynamical balance between turbulent and magnetic 
stresses within this convective layer corresponds to flux densities in excess 
of 10 15 G. Material accreting onto the core is heated by neutrinos and also 
becomes strongly convective. We compare the expected magnetic stresses in 
this convective 'gain layer' with those deep inside the neutron core. 

Buoyant motions of magnetized fluid are greatly aided by the intense 



in 
o 

o 
6 

neutrino flux. We calculate the transport rate through a medium containing 
free neutrons, protons, and electrons, in the limiting cases of degenerate or 
non-degenerate nucleons. Fields stronger than ~ 10 13 G are able to rise 
through the outer degenerate layers of the neutron core during the last stages 
of Kelvin- Helmholtz cooling (up to 10 seconds post-collapse), even though these 
layers have become stable to convection. We also find the equilibrium shape of 
a thin magnetic flux rope in the dense hydrostatic atmosphere of the neutron 
star, along with the critical separation of the footpoints above which the rope 
undergoes unlimited expansion against gravity. The implications of these results 
for pulsar magnetism are summarized, and applied to the case of late fallback 
over the first 10 3 — 10 4 s of the life of a neutron star. 
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1. Introduction 

Young neutron stars are believed to be strongly magnetized, based on measurements 
of their spin behavior ( [Kulkarni 1992| , and references therein) and the polarization of their 



pulsed emissions ( [Lyne fc Manchester 198S ). The inferred dipole fields are in the range 
10 11 — 10 14 G. In a broader sense, neutron star magnetic fields of this magnitude are weak. 
Consider, for example, a newly formed neutron star within which neutrinos of all types are 
temporarily trapped: near the neutrino photosphere (neutrinosphere) this dipole field makes 
only a tiny contribution to the hydrostatic stresses, -E>J ipole /87rP ~ 10~ 9 (£>di P oie/10 12 G) 2 , as 
compared with ~ 10~ 6 in the Sun. Over the last several years there has been a growing 
accumulation of evidence pointing to the existence of neutron stars in which the magnetic 
fields are so strong as to be transported rapidly throughout the core and deep crust of the 
star, even at the remarkably short age of ~ 10 4 yr. This evidence has come in several forms: 
the short but dramatically bright outbursts of the Soft Gamma Repeaters; the persistent, 
pulsed output of X-rays from these sources and the non-bursting Anomalous X-ray pulsars; 
and the rapid braking of the rotation of the SGRs and AXPs. These developments are 
exciting, in part, because direct physical measures of strong magnetic fields in isolated 
neutron stars such as radio pulsars have been difficult to obtain. 

The physical origin of magnetism in neutron stars is difficult to pin down because there 
are few direct observational probes of the evolving cores of massive stars. Nonetheless, a 
simple physical argument points directly to the violent convective motions that develop in 
the supernova core as the predominant source of free energy. The kinetic energy carried by 
this convection (both inside the neutrinosphere of the nascent neutron star, and between the 
neutrinosphere and the shock) is ~ 100 times greater relative to the gravitational binding 
energy than it is during any convective episode driven by nuclear burning ( Thompson 



Duncan 1993| , hereafter TD93). From this perspective, the strength of the seed magnetic 



field (present, e.g., at the end of the main sequence evolution) becomes of secondary 
consequence to the final disposition of magnetic fields in the neutron star. There is growing 
empirical evidence - from chromospheric and coronal activity in fully convective stars 
(Delfosse et al. 1998; [Kiiker fc Rudiger 1999|) and from the small scale, intranetwork 



magnetic field of the Sun (e.g. Purney, de Young, fc Roxburgh 19931) ~~ that magnetic and 



turbulent stresses reach a dynamical balance in a fluid of very high magnetic Reynolds 
number, even if that fluid is slowly rotating. 

In this paper, we consider afresh the transport of magnetic fields in the collapsed core 
of a supernova. We focus in particular on the role of the intense neutrino flux in facilitating 
transport of magnetic fields across convectively stable layers that would otherwise bury 
dynamo-generated magnetic fields. Before the supernova shock escapes to large radius, the 
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neutrinosphere is enveloped by a fairly thick layer of material (AM ~ 0.01-0. 1M Q ) which is 
stabilized to convection by a gradient in electron fraction Y e (Keil, Janka & Mueller 1996, 
hereafter |KJ1VI| ). Even as the shock propagates outward, a reverse shock forms that can 



deposit up to ~ 0.1 M Q onto the neutron core over ~ 10 3 s ( |Woosley fc Weaver 1995| ; Fryer, 
Colgate, & Pinto 1999, and references therein). 

The plan of this paper is as follows. Section 2 and the Appendix give a fresh analysis 
of persistent convection inside the neutrinosphere of a newly formed neutron star. We 
focus on the forcing of convection by secular neutrino cooling, when the lepton number 
gradient begins to stabilize the star but a negative entropy gradient persists. We connect 
our calculation to those of Lattimer & Mazurek (1981) and Reisenegger & Goldreich (1992), 
who calculate the Brunt- Vaisala frequency in hot and cold neutron stars. We also note 
that continuing fallback onto the neutron core will induce rapid convection outside the 
neutrinosphere that can deposit a strong (but tangled) magnetic field in the surface layers 
of the star. 

The rate of radial transport of a magnetic field in and around the forming neutron 
core is calculated in detail in Section |3|. Turbulent pumping moves magnetic flux outward 
rapidly through the convective core, in the direction of decreasing turbulent diffusivity 
(§ [3~1~D Outside the neutrinosphere, neutrino absorption increases the entropy of a magnetic 
flux rope with respect to its surroundings and allows it to move outward through a stably 
stratified atmosphere (§ |3.2|) . At late times, when the optically thick interior of the star 
is no longer convective, neutrinos induce transformations between protons and neutrons 
that allow a strongly magnetized parcel of fluid to move across a stablizing lepton number 
gradient (§ |3.3| ). The net result is that magnetic fields stronger than 10 13 — 10 14 G will rise 
buoyantly to the surface of a newly formed neutron star. The equilibrium configuration of 
a magnetic flux rope that extends outward from the convective interior of the star through 
the dense outer atmosphere is addressed in Section 4. 

We summarize in §5 the implications for pulsar magnetism, discussing the 
relative importance of the interior and exterior dynamos, and applying our results 
to the case of late fallback (over ~ 10 3 — 10 4 s). The gravitational binding energy 
of the accreting material is converted to electron-type neutrinos, whose luminosity 
GMM NS /R NS = 7 x 10 49 (M/lO" 4 M s- 1 ) (M NS /1.4M ) (.Rns/10 km)" 1 erg s" 1 is smaller, 
by about two orders of magnitude, than the luminosity emitted during the prompt ~ 30 
s Kelvin-Helmholtz phase. Nonetheless, this neutrino flux is high enough to induce many 
convective overturns within ~ 100 km of the neutron star (Thompson 2000, hereafter [I'OCj ; 
§H), and thus to wind up a seed magnetic field within the accretion flow. As the result of 
the 
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2. Amplification of Magnetic Fields 

The supernova core is subject to a violent convective instability inside the v-sphere, 
driven by optically thick neutrino cooling; and a distinct convective instability below the 
shock, triggered by heating by the outward-streaming neutrinos (Fig. 1). This second 
instability will persist as long as substantial accretion continues onto the neutron core, and 
the core remains a luminous source of neutrinos. We first show that entropy gradients will 
drive convection below the neutrinosphere of a young neutron star. We then consider the 
response of a magnetic field to instabilities either above or below the neutrinosphere. 



2.1. Convection Inside the Neutrinosphere 

The stability of a neutron star to convection depends on gradients in both entropy 
and composition. The key compositional parameter is the net lepton number per baryon, 
Yi = Y e + Yy e , where Y e and Y Ve are respectively the numbers of electrons and electron 
neutrinos per baryon. During the first ~ 10 seconds, when the star is optically thick 
to neutrinos, it is unstable to convection ( |Wilson fc Mayle 1988] ; Burrows 1987; TD93; 
K.IM: [Pons et al. 1999| ). The Ledoux criterion for convective instability is expressed most 



conveniently as 

dS (dp/dY^g dYt _ dS 
dR (dp/dS) PXl dR dR 

(e.g. |Lattimer fc Mazurek 19811) . 




It is worth explaining here why this deep convective instability is generic, because its 
existence has been called into question on the basis of hydrodynamical simulations of the 
outer neutron core, at relatively low optical depth to neutrinos flMezzacappa et al. 1998[ ). 



Deep convection is distinct from the short-lived overturn which occurs where the weakening 
shock establishes a negative entropy gradient and which, according to some simulations such 
as those of [Burrows, Hayes, fc Fryxell (1995)| and Bruenn & Mezzacappa (1994), may be 



quenched by neutrino diffusion in the first 20-50 ms after bounce. In deep convection, the 
gradients of Yi and S both tend to be reduced (e.g. |KJM]) , but are not entirely eliminated. 



Rather, they are constantly regenerated by radiative losses from the neutrinosphere. Thus, 
convection must continue until the magnitude of the radiative gradient drops below the 
adiabatic gradient. The reduction of the radiative gradient results primarily from the 
decrease in flux over the Kelvin time, which is of order ten seconds. 



Ledoux convection is always encouraged by negative dS/dR, but the effect of dYi/dR 
on convective instability depends on the temperature and composition. Negative dY\jdR 
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is quickly established in a cooling neutron core by neutrino transport at high Y\ and by 
chemical equilibrium between electrons, protons and neutrons at at low ^ 0.1 (where the 
lepton number is carried almost entirely by the electrons). The sign of the thermodynamic 
derivative (dS/dYi)p yP is therefore crucial to the nature of the convective instability. The 
negative lepton-fraction gradient is de-stabilizing [(dS/dYe)p, p < 0] if the neutrino chemical 
potential fi Ue is high enough and if the temperature is lower than a critical value T*. In 
this situation, where a compositional gradient is potentially the main driving force behind 
the convective instability, salt-finger effects may play an important role in its non-linear 
development fl Wilson fc Mayle 19881) . 

The critical temperature T* can be estimated analytically when the protons, neutrons, 
electrons, and electron neutrinos are all treated as ideal and almost degenerate fermi gases. 
As derived in Appendix A. 2 (eq. | |A35| j), one hasQ 



35 



He 




-1/2 


3(2fi Ve - He) 


1/2 


7T 






m n c 2 





1/2 



MeV (Y e ~ 0.1) 



(2) 



Here fi e and fi Ue are the chemical potentials of the electrons and electron neutrinos, 
respectively. We have scaled the density to the nuclear saturation density n sat = 1.6 x 10 38 
cm -3 (corresponding to a mass density m n n sat = 2.8 x 10 14 g cm -3 ). This critical 
temperature is comparable to the initial peak temperature in the neutron core. Negative 
dYi/dR indu 
equivalently 



dYi/dR induces convection below the temperature T*, but only as long as fj, Ue > |// e , or 



Y„ 



>ly, 

16 



(3) 



In fact, Y Vr falls below this bound quickly in the outer parts of the neutron core, in about 1 
second according to the cooling models of Pons et al. (1999). Thereafter, the composition 
gradient becomes stabilizing ( [Lattimer fc Mazurek 198~H |Goldreich fc Reisenegger 1992j 
KJM| ; see §A.4| ). However, we now show that a negative entropy gradient can drive 
convection in a multicomponent, noninteracting Fermi gas, even in the presence of a 
stabilizing composition gradient. 

Such a "late phase" of convection occurs if the radiative temperature gradient is 
sufficiently super-adiabatic. The radiative gradient is set by the neutrino opacity, the energy 



1 Neglecting terms of the order of Y Ve /Y e and /i e /m„c 2 , but not fj, Ue /fj, e = (2Y l , e /Y e ) 1 ^ 3 . 
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flux, and the temperature: 

(7K\ 16as B T 4 / 1 \ dWT_ 
rad " v 8 7 3 ^Kff^j dlnR- 1 J 

Here, cr„ is the neutrino cross section and the angular brackets indicate a Rosseland mean. 
During the intermediate stages of cooling, the energy flux is carried primarily by \l and r 
neutrinos: their opacity is due to scattering and is smaller than the absorption opacity of v e 
on neutrons above a temperature of ~ 10 MeV (where the absorption is not Fermi-blocked 
by the degenerate electrons; Pons et al. 1999| ). In this situation, the effective number of 



neutrino species in expression (U) is N v ~ 2. However, the production rate of fi and r 
neutrinos may be so low that they cannot carry the flux, in which case N v = 1. 

A negative entropy gradient is generated by steady radiative transport, as is easily 
seen from the following argument. The opacity of non-degenerate neutrinos scattering from 
degenerate nucleons scales as E% (as compared with the E 2 . energy dependence of the weak 
neutral current in a non-degenerate plasma; [Iwamoto 1981[ ). The Rosseland mean free path 



is (cr)(T) ~ <j(E u = 2ksT) for a Fermi-Dirac distribution of neutrinos with temperature T 
and vanishing chemical potential. The mean-free path for neutron scattering is 

1 —— = 2 x 10 3 ( n E " \ cm (5) 

n n a(E v ) V30 MeV/ v ; 



at a density n b = n sat (we normalize to Fig. 12 of [Reddy, Prakash, fc Lattimer 199"8"| ). The 



radiative temperature gradient can, then, be directly related to the total radiative flux 
F rad = Lv/AirR 2 : 

,diaT " -i.3 m i (^^) r^r (^r. o 



dlnRJ iad \2J \1Q 52 erg s" 1 / V10 MeV J V10 km 



This radiative gradient must be compared with the gradient required to trigger 
convection in the face of a (possibly) stablizing composition gradient. This convective 
temperature gradient is calculated by setting dS/dR = (dS/dYi) PtP dYi/dR. The entropy 
(per baryon) is provided mainly by free neutrons and protons, with respective abundances 



Y n ~ 1 - Y e and Y p = Y e : 



k 



S n i , . - i " / ,-■ i /^x 



B 



Y e ( —] + (1 " Y.) | ^ 



Here, the fx p , \i n are the chemical potentials of the degenerate protons and neutrons 
respectively. This critical temperature gradient for convection is calculated in Appendix 
A 3- 

fdlnT\ 2 fdlnp\ Y^ 
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m n (kTy 




d\nY e 
d\nR' 



Notice that the second term (proportional to din Y e /d In R) is guaranteed to be stabilizing 
when din Y e /d\nR < and 2fi Ue < \i e . The coefficient 

m n (A;T) 2 i ' U U sat Ao.l J ^20 MeV J ' W 

This analytical criterion indicates a slightly earlier end to convection than is found 
in the cooling models of Pons et al. (1999) (cf. their Figs. 9, 21, and 24). The value of 
(d In T/d In r) conv depends on a cancellation between two terms of opposing signs (on the 
RHS of eq. ||). At a density rib — n sat and a density gradient d In p/d In r ~ —3, the 
radiative temperature gradient becomes too weak to induce convection below a luminosity 
L y ~ 1 x 10 52 erg s _1 and temperature of ~ 20 MeV. Since the binding energy is 3 x 10 53 
ergs, the luminosity will be high enough to drive convection for 10 seconds or more. 

Unlike a main sequence star, the thermal energy in the young neutron star is not 
generated in the central core. Heat is generated through the entire star as the lepton 
fraction relaxes to its equilibrium value. The absence of shock heating in the inner half of 
the star leads to a positive temperature gradient at small radii, which delays the onset of 
convection deep in the star ( [Burrows fc Lattimer 1986| ). However, as the outer layers cool 



by neutrino loss, a negative entropy gradient is produced at smaller and smaller radii. 

The above argument shows that passive neutrino transport through a newly formed 
neutron star leads inevitably to convective instability. Once the inconsistency of purely 
radiative heat transport is established, one knows that some fraction of the heat flux must 
be carried by convection. In other words, the negative entropy gradient is not entirely erased 
by the convection (as is apparent in the simulations of KJM). Rather, it is maintained at a 
large enough magnitude to carry the flux needed to supply that lost at the neutrinosphere 
- until the magnitude of the radiative gradient drops due to decreasing luminosity or 
increasing neutrino mean free path. 

The convective energy flux is easily estimated using mixing length theory. In terms of 
the convective velocity V con , one has 

i^con pKon- (10) 

Numerically, V con ~ 10 8 cms _1 , which is much less than than the speed of light and 
substantially smaller than the sound speed c s ~ 5 x 10 9 cms _1 . Notice that the convective 
velocity varies as and so depends weakly on the convective efficiency. As we review in 
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the next section, this means that the maximum possible dynamo- generated magnetic field 
will also vary weakly with the convective efficiency. 

The latest numerical simulations reach differing conclusions about the persistence of 
Ledoux convection inside the ^-sphere, but for reasons that can be ascribed in part to the 
handling of neutrino transport. KJM employ a fully two-dimensional hydrodynamic code, 
with neutrino transport restricted to radial rays, and find that convection continues over the 
entire Kelvin time. The simulations of Mezzacappa et al. (1998) further restrict neutrino 
transport to one spatial dimension, and find that the two-dimensional convective motions 
die out after a few tens of milliseconds. They argue that Ledoux convection is inhibited 
by neutrino transport in the outer layers of the neutron core (p <^ 10 13 g cm -3 ), where the 
diffusion time ~ T Ue £p/c is shorter than the convective overturn time in the absence of local 
neutrino heating and cooling. However, deeper in the core the convection overturn time 
Tcon = £p/V con ~ 1 — 3 ms is much shorter than the Kelvin time and this conclusion is 
reversed. The cooling models of Pons et al. (1999), which include more realistic neutrino 
opacities, find that a deep convective instability persists for tens of seconds. 

A recent paper by Miralles, Pons, and Urpin (2000) carefully treats how convective 
instability is modified by the transport of heat and lepton number, and by viscosity. They 
find that the cooling neutron star simulated by Pons et al. (1999) is subject to a slower 
double-diffusive instability (the 'neutron finger instability'; [Wilson fc Mayle 1988| ) even 
in regions that are stable according to the Ledoux criterion (|l|). Nonetheless, at both 1 s 
and 20 s, the region of the star that is subject to an unstable g-mode with a fast growth 
rate 1 msec; Fig. 3 of Miralles et al. 2000) closely approximates the region found 
previously to be unstable according to the Ledoux criterion (Fig. 24 of [Pons et al. 1999[ ) . 



The new 'convective' instability criterion (eqs. [25] and [26] of Miralles et al. 2000) is more 
restrictive than Ledoux, so that regions of the star that are unstable according to eq. ([!]) 
are labelled as being subject only to the 'neutron finger' instability (in spite of the equally 
fast growth rate). In fact, one can show that this new 'convective' instability criterion 
remains more restrictive even in the limit of slow transport, where the Ledoux criterion 
should be accurate. For that reason, we interpret the results of Miralles et al. (2000) to 
imply that transport effects do not significantly restrict the region of a newborn neutron 
star that undergoes a vigorous convective instability. 



2.2. Dynamical Balance between Convective and Magnetic Stresses 



Even in the absence of rotation, these convective motions carry enough energy to 
amplify a magnetic field to enormous strengths. The argument for such a stochastic dynamo 
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is based on a phenomenological scaling, rather than purely theoretical considerations, and 
is presented in detail in TD93. One observes in the convective envelope of the Sun a small 
scale ('intranetwork') component of the magnetic field, with a mean pressure about 10 
percent of the turbulent pressure, 

= e B -0.1. (11) 

r v con 

This field is present even at Solar minimum, when it is distributed almost uniformly 
across the Solar disk (|Murray 1992|). It has been conjectured to be the consequence of a 



stochastic, non-helical dynamo that builds up magnetic flux on the scale of an individual 
convective cell (TD93; purney, de Young, fc Roxburgh 1993| ). Such a dynamo would 



operate almost independently of the global Solar dynamo that manifests itself at the surface 
in the form of sunspot activity and the dipolar magnetic fieldQ. Some theoretical support 
for this conjecture comes from models of dynamo action in mirror symmetric turbulence 
( Ruzmaikin fc Sokolotf 1981] ; |Meneg uzzi fc Pouquet 1989| ). 



The total convective kinetic energy grows from E con — 0.5 to 1 x 10 49 ergs between 
0.1 and 1 s, as the convective zone expands to encompass the entire neutron star (KJM). 
At 0.1-0.2 s after bounce, a zone interior to 0.9M Q is unstable to Ledoux convection, with 
overshoot extending out to an enclosed mass of 1.05M Q . This is in accord with a simple 
mixing length argument (Purrows 19871 TD93). The convective luminosity is related to the 
luminosity L u = L u52 x 10 52 erg s _1 in all neutrino species via L con — L u — AnR^ on pV^ on . 
The kinetic energy in convection E con = -f-R^ x \pV^ on can then be expressed in terms of 
the convective overturn time 

Tcon = i P /V con = 2 X 10- 3 ( — r ) S (12) 

V30 km/ VSxlO^^cms^ 1 / 

via 

E cori = Wcon ~ (\ ~ l) X 10 5 ° (^) (^-) ergS. (13) 

The pressure scale height lp — R/5 while the nucleons are only mildly degenerate. Note 
that the convective overturn time is short compared with the expected initial spin periods 
of most neutron stars ( |Duncan fc Thompson 1992[ ). Applying the scaling ( TTTT ) for the 



2 The two dynamos are not, of course, entirely independent: the global dynamo feeds magnetic flux into 
the convection zone. However, one infers from the time- variation of Solar p-mode frequencies with Solar cycle 
that the proportionality between magnetic and turbulent stresses is maintained down to least six pressure 



scaleheights below the photosphere ( Goldrcich ct al. 1991 ). This strongly indicates a local amplification 



mechanism that leads to a dynamical balance between the two stresses. 
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magnetic energy implies a total magnetic energy 

s > s rs d, '= !A - u "©(i^) < i4) 

inside the v-sphere. The corresponding r.m.s. magnetic field is 

2)T = „ 15 x 10 „ (_^) 1/2 f^y m G . (15) 

V lRv V R 3 J VO.iy V 10 50 ergs/ V30W v ; 

It should be noted that this dynamo process taps the energies of the diffusing 
and u T , as well as the is e , v e . By contrast, direct neutrino heating outside the //-sphere is 
mediated by absorption of v e and v e only 

The region below the shock provides a second potential site for magnetic field 
amplification. Heating by the outward flux of electron-type neutrinos, through the 
charged- current absorption u e (u e ) + n{p) — > e~(e + ) +p(n), induces strong convection 
(e.g. Burrows et al. 1995; |Janka fc Mullcr 1996|) . The conditions are especially favorable 
for magnetic field amplification if prolonged infall continues after the shock succeeds 
(leading to the formation of a secondary accretion shock) or if the expansion of the 
shock is asymmetric, so that accretion continues in one hemisphere even as the shock 
expands in the opposing hemisphere ( |1'00| ). Buoyancy forces will pin magnetized material 
just below the shock even while heaver material settles below the shock toward the 
neutrinosphere. The convective Mach number is expressed in terms of the accretion rate 
M, shock radius _R s h, core mass M core , and the (normalized) electron neutrino luminosity 
C = Y n L De «4 e >/200 MeV 2 ) + Y p L Ue («>/200 MeV 2 ), 

(M\ n« V3 (K\ 1/2 ( R sh y 1/6 /£ 52 \ V3 / Mcore V 1/2 /-i K \ 

(MaU = 0.6 e h / at y (^^j [—) (^ TJ ^-j . (16) 

Here £heat is the dimensionless ratio of the net heating rate (from which cooling 
by electron and positron captures has been subtracted) to the gross heating 
rate (uncompensated by neutrino cooling). The convective overturn time is 
r con ~ 2 x 10~ 3 £i^ 3 (iU/100 km) 5 / 3 (C 52 /A)-^ 3 s (§2.3 of p0|) . This overturn is 
fast enough to allow many efoldings of the buoyant magnetic field, if the shock radius -R sn 
has collapsed to ~ 100 km. When the magnetic stresses reach a dynamical balance with 
the convective stresses, according to equation (|TT|), the strength of the field has increased 
to ( [TOOD 

{B} " 1Xl ° £b -i £ ^ [m q s-i) UooW U ) {l0j \1AM Q ) 

(17) 
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2.3. Convective Dynamo During Late Infall 

It is also worthwhile estimating how this equipartition field will diminish, as the 
rate of accretion drops off. Late infall (AM = O(10 _1 ) M over ~ 10 3 — 10 4 s) is a 
plausible consequence of the formation of a rarefaction wave in the exploding supernova 
core ([Woosley fc Weaver 1995| ; Fryer, Colgate, & Pinto 1999). Let us suppose that M 
remains large enough that neutrino cooling at the base of the accretion flow continues to be 
dominated by electron/positron captures on free protons and neutrons, as compared with 
annihilation e + + e~ — > v + v (as assumed in the models of Chevalier 1989). We focus 
on the inner regions of the flow, R ~ 30 — 100 km, well inside the accretion shock. Here the 
temperature is high enough for helium to be photo-dissociated, but the capture rates on free 
n, p are still small enough to prevent the formation of a stabilizing gradient dY e /dR. The 
thermal pressure is then dominated by photons and (at least mildly) relativistic pairs, 
P = ^faT 4 . The corresponding isentropic hydrostatic profile is T(R) oc R~ x , p(R) oc -R -3 , 
and P{R) ~ \{GM NS / R)p{R) oc iT 4 . The base of this settling flow may be defined as 
the radius where the flow time (\R/Vr)r across a pressure scale-height equals the time 
(3P/Q eJV )i? for captures to remove the thermal energy. (The neutrino emissivity per 
unit volume scales as Q eN oc T 6 p in a medium with non-degenerate electrons.) Cooling is 
slower than advection outside this radius, and the gradient dY e /dR is too weak to stabilize 
the flow against the effects of neutrino heating ( [Mezzacappa et al. 1998| ). As a result the 
flow becomes violently convective. 

Combining the constraint of mass conservation, M = 4tt R 2 VrP, with the equation of 
hydrostatic balance at the base of the settling flow, and taking the cooling radius to be a 
fixed multiple of the stellar radiusQ, one finds the scaling relations 

p(R) oc M 2/5 ; T(R) oc M 1/w (18) 

at a fixed radius R. 

From the relations (Jl8|), one can deduce several things: 

1. The temperature at the base of the settling flow typically exceeds m e c 2 , for 
parameters of interest, as is needed for the pressure to be dominated by photons and 
relativistic e ± pairs. For example, the temperature is ~ 1 MeV at the base of an accretion 
flow carrying a mass flux M ~ 10~ 4 M Q s -1 . 

2. The nucleons in the settling flow become less degenerate with decreasing accretion 



3 In the presence of an extended atmosphere, the stellar radius is conveniently defined as the radius where 

P = Psat ■ 
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rate: the ratio of temperature to fermi energy scales as kgT ' je v oc M" 1 / 6 . At the same 
time, the thermal energy per baryon remains almost constant at fixed radius, T A j p oc M°, 
and so effective dissociation of helium is expected at the base of the accretion flow, even at 
relatively low accretion rates. 

3. At an intermediate infall rate ~ 10~ 4 — 10 -5 M s _1 , the assumption of neutrino 
cooling by captures is indeed self-consistent (as compared with e ± annihilation 

e + + e~ — > v + v). The relative rates of energy loss through these two channels vary weakly 
with M at a fixed radius. In particular, Q eN /Q e± oc p/T 3 oc M 1//10 , before e ± captures have 
removed enough heat to make the electrons degenerate. 

4. Neutrino heating still causes the infalling material to be violently convective, 
before finally settling onto the neutron star surface ( TOO ). At late times, the luminosity in 
electron-type neutrinos can be directly related to M, 

GMM m __ in49 / M \ ( M m \ ( R m y 1 _! 



^ + L °> = ^ = 7 x 1048 (TOW J ) iToliJ ■ ' 19 > 

The heating rate (per unit volume) due to charged-current absorption of these neutrinos 
on the settling baryons is 3P/rh ea t oc pM (eq. [P7[]). The energy density in the resulting 
convective motions scales as \pV^ on ~ (3P/r heat ) x (£p/V con ). As a result, the convective 
speed decreases more slowly with decreasing accretion rate, V con (R) oc t^J 3 oc M 1//3 , than 
does the settling speed, Vr(R) oc M 3//5 . Moving outward through the accretion flow, the 
capture (cooling) timescale grows as ~ p _1 T -2 ~ R 5 , given an approximately isentropic 
hydrostatic profile T oc and p oc R~ 3 . By contrast, the heating time increases only as 
r hcat ~ R 2 - Thus, on general grounds one expects a 'gain region', within which dS/dR is 
strong enough, and dY e /dR is weak enough, to trigger convection (e.g. Burrows et al. 1995; 
Janka fc Muller 1996|) . 



5. Finally, using the above scaling relations, it is easy to work out how the equipartition 
magnetic field scales with accretion rate: 

(BY 2 = z¥(^P) 1/2 V CO nOcM 8 / 15 , (20) 

at a fixed radius. The equipartition field drops by a factor ~ 0.01 from M = 1 M s _1 to 
M = 1O" 4 M s _1 . 



3. Transport of Magnetic Fields 



We now discuss the transport of a magnetic field through the convective interior 
and convectively stable exterior of a nascent neutron star. The star's convective core is 
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surrounded by an overshoot layer of thickness AM ~ 0.1M© and another convectively stable 
layer of similar mass that extends out to the gain radius (KJM). Magnetic fields amplified 
in the convective region are transported rapidly into an adjoining overshoot layer (with a 
velocity ~ V con /£p). Recent simulations of magetoconvection ( |l'obias et al. 199"S|) show the 
the magnetic field will remain pinned at the base of a convective layer, even against the 
action of buoyancy forces. Thus, we expect the overshoot layer above the convective core of 
a neutron star to acquire magnetic fields as strong as those in the core. Transport of the 
magnetic field through the surrounding convectively stable layer is greatly aided by i) the 
concentration of the field into a dense fibril state and ii) rapid neutrino heating. We also 
discuss an additional subtlety involving turbulent pumping, driven by a radial gradient in 
the turbulent diffusivity. 



3.1. Turbulent Pumping within Convective Zones 

A nascent neutron star and a late-type main sequence star both are deeply convective. 
In the outermost convective layers of a main sequence star, the turbulent diffusivity 
increases outward due to the recombination of hydrogen, which buffers the decrease of 
temperature with decreasing density. By contrast, there are no similar ionization effects 
that enhance the neutrino scattering or absorption cross sections in a neutron star and, as 
we now show, the turbulent diffusivity 

X = \v co J P (21) 

decreases outward. Once the ^-sphere shrinks to ~ 15 — 20 km, the pressure scale 
height with the surface layer of mildly degenerate nucleons scales with density as 
tp = P/pg oc p 2 / 3 R 2 /M(< R), where M(< R) is the mass enclosed within radius R. 
Constancy of the convective energy flux implies that V con oc p _1 / 3 i?~ 2 / 3 . Combining these 
two expressions yields 

1/3^4/3 

* * W (22) 

One sees that \ decreases with radius in the outer part of the convection zone, where p 
decreases at almost constant R and M(< R). A similar conclusion holds for the convective 
shell in between the gain radius and the shock. 

A passive contaminant such as smoke, or a magnetic field, experiences a net drift 
velocity in the direction of decreasing turbulent diffusivity, Vdrift = — V(V con £ p ), whose 
magnitude is comparable to the convective velocity, | Vdrift | ~ V con . Magnetic flux ropes are 
convected downward in the upper layers of the Solar convection zone, and their buoyant 
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rise to the surface is prevented, as long as thrift ^ (a/^p) 1 / 2 B/^/Anp (e.g. Ruzmaikin & 
Vainshtein 1978). (Here 2a is the width of the flux rope.) By contrast, the mean drift 
velocity is directed upward at the top of a neutron star convection zone, and the magnetic 
heating rate should be correspondingly higher. We emphasize that this effect is more 
pronounced if the flux ropes are tied together into a larger network, so that the stochastic 
motions of individual flux ropes within individual convective cells are suppressed with 
respect to the mean drift velocity. 

A related effect, which has been studied in the context of Solar magnetism, involves the 
asymmetry between cold downdrafts and hot updrafts in a compressible, convective layer. In 
the Solar convection zone the downdrafts are much denser and narrower than the updrafts 
and, as a result, have a higher vorticity. Radiation- hydrodynamical simulations ( porch| 
fe Nordlund 2000| ) show that the magnetic field lines are wound up in the downdrafts, 
and suggest a tendency for the field to be pumped toward the bottom of the convective 
layer. The convective region of a newly formed neutron star is, however, supported by the 
pressure of mildly degenerate neutrons, and is less compressible than the Solar atmosphere. 
Numerical simulations (KJM) suggest a greater symmetry between downdrafts and updrafts 
than is inferred for the Sun. 



3.2. Transport Aided by a High Neutrino Flux 

In this section, we examine how the buoyant transport of a magnetic field through a 
supernova core is facilitated by the intense neutrino flux out of the core. The neutrinos 
have two principal effects, which we examine in turn. The first involves heating of the 
non-degenerate material outside the v-sphere, which equilibrates the temperature between 
a flux tube and its surroundings, and allows the tube to rise buoyantly across a convectively 
stable layer. The second effect involves transport of a magnetic field through the degenerate 
neutron core after convection has stopped and the core has become stably stratified. In such 
a situation, the transport rate is determined by the speed with which the electron fraction 
returns to its beta-equilibrium value QPethick 1992| ; |Goldreich fc Reisenegger 1992|) . Long 



after the neutron star has radiated most of its original heat, the rate of modified URCA 
reactions may be kept high by magnetic dissipation (Thompson and Duncan 1996). During 
the first few seconds of the life of a neutron star, direct URCA reactions are possible, and 
proceed at an extremely high rate dominated by charged current absorption of v e and v e on 
free nucleons. 



We now focus on the early stages of a supernova explosion, before the success of the 
shock, and the onset of a low density, neutrino-driven wind outside the v-sphere. The 



- 15 - 



calculation separates into two cases, depending on whether the pressure is dominated by 
nucleons or by radiation and relativistic pairs. Curvature and tension forces, discussed in §4, 
are now ignored and the rope is assumed to sit horizontally in a plane-parallel atmosphere. 



3.2.1. A Nucleon Pressure-Dominated Zone Outside the Neutrino sphere 

The outer atmosphere of a newly formed neutron star is convectively stable, with 
entropy and electron fraction both having positive radial gradients, dS/dR, dY e /dR > 0. 
Free neutrons within this layer absorb electron neutrinos at the rate 

E v +n ^ p+e - = ^ ( ^) ( {£ :^ r2 ) ( — 1—) 2 MeV neutron- 1 s" 1 (23) 



f V 4 J \200 MeV 2 / V100 km 

( [Bethe fc Wilson 1985| ). The dependence on the mean square neutrino energy is derived 
from the weak interaction cross section. The expression for proton heating (p e +p — > n + e + ) 
is obtained by replacing L Ue and {e 2 Ue ) with L Pe and (e 2 ^). Here / ~ | for isotropic emission 
at the ivsphere, increasing to / = 1 further out where the neutrinos free stream. The net 
heating rate due to absorption of v e j v e on free neutrons (Y n per nucleon) and protons (Y p 
per nucleon) is 

Q = (^t-^'^e+n^p+e- 4" YpEp e +p— >n+e~ ) j (24) 

per unit volume. Heating is compensated by neutrino cooling driven (primarily) by electron 
and positron captures, so that the net heating rate becomes Q uN + Q eN = Q vN [\ — (T/T eq ) 6 ]. 
The equilibrium temperature is 

Here, we define the convenient dimensionless parameter 

£52 = Yn Lv ° 52 (2OO MeV 2 ) + Yp Lpe 52 (2OO MeV 2 ) ' (26) 



The time for the temperature to relax to the equilibrium (p5|) is r heat ~ (^pk B T/m n )/6Q 



2' 

or 

r heat -bxlU f{ A ){ 5Qkm ) { 4 ) [ 2QQ MeV2 ) 



(27) 



since the heat capacity is dominated by free nucleons. The factor of | arises from the T 6 
dependence of the e ± -capture cooling rate. The effectiveness of neutrino heating can be 
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measured by comparing r heat with the dynamical time r dyn = £p/c s ~ (R/5)(5k B T /3m n ) x l 2 : 

Theat 1 Rfrr S/2 ( R \ ( Y n L Ue 52 + Y^Lp e 52 \ 1 / ) \ , , 

^- L6/T4 loOkmV I 4 J 1,200 MeV^ J " (28) 

This implies Th Ga t ^ 7d yn at the i/-sphere (i?^ ~ 50 km at 0.1 ms, decreasing to R u ~ 20 km 
at 1 s; e.g. |Mayle fc Wilson 19891 ). 

The nucleon pressure is suppressed inside a magnetic flux rope that is in pressure 
equilibrium with its surroundings (Fig. 2.), 

B 2 pMk B T M _ pk B T . . 

- — i . [ZU 

87r m n m n 

However, the pressure of the electrons is only weakly perturbed by the magnetic field 
when Y e is small. The density is high enough that the electrons are moderately degenerate 

Oe/T £ 7r), 

P » ^ k fflv n = 2 x 108 Y v l T McV g cm" 3 . (30) 

At a fixed radius, the equilibrium electron pressure is the result of a competition between 
absorption of u e /p ei 

^ ex R~ 2 [L Ue (e Ue )n n - L De {e, e )n p ] ~ R~ 2 L„ e {e„ B ) (1 - 2Y e ) , (31) 

and the capture of (degenerate) electrons on nucleons, 

^ oc -£p oc -Y 2 p\ (32) 

The electron fraction therefore varies with density as 

1 - 2Y 



Y 2 



oc p 2 , (33) 



at fixed R. The decrease AP e in electron pressure accompanying a decrease Ap in density 
is small where F e <l; 

Pe 3(1 -n) p 1 ' 

Furthermore, P e is only a fraction of the nucleon pressure because Y e is small where p is 
large: 

Pe _ ^ _ n , / g \ fTuoV\ /^\ 4/3 
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We conclude that the flux rope is less dense than its surroundings, pu < p, when its 
temperature rapidly equilibrates with its surroundings. Rapid neutrino heating, Th eat ^ T dyn , 
brings the temperature inside the flux rope close to the ambient value, T M ~ T. The 
density deficit inside the rope is then pu — p — —p/3p, where 



01 



B' 



8vmT 



(36) 



Balancing the buoyancy force (integrated over the cross-sectional area ira 2 of the rope) 
against the drag force Cd x 2aV^, yields the equilibrium vertical speed 



V R 



tt a \ V2 B 



AC, 



d^p 



(37) 



(cf. Parker 1979). 



The vertical speed Vr is determined by the rate of neutrino heating when Theat ^ T &yn- 
The force balance pu Vp/ a ~ g(p — Pm) implies that the density is almost uniform, p M ~ p, 
and that there is a temperature deficit 

T M -T = -f3 P T (38) 



inside the rope. As the rope moves upward, Tm changes due to neutrino heating (eq. p7fl ) 
and adiabatic cooling, 



1 dT, 



M 



1 T M -T, 



cq 



M 



dt 



7~heat 



T, 



M 



+ (7 " 1) 



1 dp 



M 



p M dt 



(39) 



Note that it is the equilibrium temperature (P5[) , rather than the ambient temperature T, 
that enters on the right side of (|3~9"D. Combining the conditions of pressure equilibrium and 
flux conservation, B/pm = constant, yields the following relation between time derivatives, 



1 dT, 



M 



1 dP 
~P~dE 



T M dt 

Equations (|38|), ([39|) and fllCf ) together yield 
1 (3 P T eq -(l-0p)(T-T eq ) 



(1 + 2/3, 



1 dp 



'M 



Pm dt 



Theat 



T: 



M 



F dt pm dt 



(7- 



(40) 



1) (^_2£pdp^ 



dt pm dt 



where S is the entropy per baryon in units of hp- The vertical rise time ip/Vp of the rope 
is obtained by substituting d/dt — > Vp(d/dR) and pm — P, 



V R 



7~licat 



[1-Pp)T 



(3 P T eq -(l-(3 P )(T-T ( 



cqj 



R 



(7 



dS 
dhiR 



2f3 f 



dlnp 
d\nR 



(42) 
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In the marginally failed shock model of Janka & Mueller (1996) one finds dS/dlnR ~ +12, 
d\iap/d\nR ~ —4, and £p/R ~ | within the stably stratified layer outside the v-sphere at 
100 ms after bounce. Thus, 

(l ' l-6r heat ( l -^M (43) 



v R v ft 



when T ~ T cq . 



It is straightforward to see that the upward, buoyant rise of a closed loop of magnetic 
flux will continue through the convectively stable layer outside the //-sphere and reach 
the convective shell interior to the shock. We consider both cases (3p -C 1 and (5p ~ 1 in 
turn. In the first case, the radial transport time is limited by neutrino heating (eq. | ]43| | ) 
and scales as (3p l fTR 2 oc (3p l fR. The radial dependence of (3p is obtained by noting that 
the density inside the loop scales as pu ~ p, and the flux density as ~ . Thus, 
/3p = B 2 /8nP oc p 4 / 3 /P oc i? -1 / 3 and (using the same scalings p oc i?~ 4 and P oc i? -5 ) 
the transport time increases with radius. At the inner boundary of the convective shell 
(R ~ 100 km) the transport time is only ~ 1 x lO" 3 ^ 1 {L Ue52 /A)- 1 ((4J/200 MeV 2 )^ 1 s. 

Transport is even faster when (3p ~ 1, since buoyancy is effective even in the absence 
of neutrino heating. Balancing B 2 /87r = P, and making use of the vertical velocity (]37|) of 
a horizontal flux rope, one finds 

(44 > 

This yields £p/Vp oc _R 11//8 , using the same scalings for P and p. 



3.2.2. A Radiation Pressure-Dominated Zone Outside the Neutrino sphere 

We now consider transport of a magnetic field through a convectively stable atmosphere 
whose pressure is dominated by relativistic pairs and photons. The zone just interior to 
the gain radius fits this description. The temperature of the settling flow exceeds the 
equilibrium value T eq (eq. f25|), and neutrino cooling dominates heating. The pressure 
equilibrium condition becomes 

H + l aT *' - T2 aT> - (45) 

In this situation, the nucleon density inside the flux rope is an independent variable, 
which determines the strength of the buoyancy force and the rate of neutrino heating. When 
Tm T 4 , the equilibrium magnetic flux density decreases with radius as B oc T 2 oc R~ 2 . 
The density inside a rising, horizontal flux rope decreases according to pu oc B oc R~ 2 . By 
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contrast, the density of the surrounding settling flow decreases much more rapidly — faster 
than p(R) oc R~ 3 because the relativistic entropy per baryon (S oc T 3 /p) increases with 
radius. Thus, a flux rope will traverse several pressure scale heights only if it starts off at 
a density well below that of the ambient settling flow. In such a situation, direct neutrino 
heating is not required to ensure buoyancy. In practice, an effective barrier to magnetic 
buoyancy is not expected during the first ~ 1 s, because the flow remains radiation-pressure 
dominated only within a short distance below the gain radius. 

Radiation pressure also has a significant effect on the equilibrium shape of a magnetic 
arcade with fixed endpoints. An arcade of magnetic flux, whose footpoints are separated 
by an angle A<f>, will open to infinity when A</> exceeds the critical value Q79D given below. 
This critical angle depends on the dimensionless parameter £ P oc pT/P (eq. |7Dfl), in such a 
way that the critical footpoint separation decreases as photons and pairs contribute more 
to the pressure. 



3.3. Transport in a Convectively Stable Neutron Core 

We now consider the transport of a magnetic field in the central neutron core, after 
the core has become convectively stable. The electrically charged particles (protons and 
electrons) are tightly coupled to the magnetic field on short timescales. Transport of these 
charged particles through the degenerate neutron fluid proceeds extremely slowly, and at 
high temperatures is limited by neutron-proton drag (Goldreich and Reisenegger 1992). As 
a result, the interior of a newborn neutron star is very nearly an ideal magnetofluid. For 
this reason, we consider the bulk hydrodynamical motion of an isolated rope of magnetic 
flux. As in the previous section, we neglect curvature and tension forces and assume plane 
parallel symmetry. 

The neutron core becomes convectively stable as it cools because /^-equilibrium value 
of Y e has a negative gradient, and the density of cool degenerate matter decreases with 
decreasing Y e ( [Reisenegger fc Goldreich 1992| ; § [A.4| ). This happens, at the latest, when 



the temperature drops sufficiently that neutrinos begin to escape directly from the core. 
A flux rope that starts out in /^-equilibrium is less dense than the surrounding fluid; as 
it rises buoyantly, it falls out of /^-equilibrium and its density equilibrates. We treat the 
neutrons, protons and electrons as three ideal, degenerate Fermi fluids. The departure from 
/3-equilibrium is measured by the chemical potential imbalance 
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In pressure equilibrium, 



B 2 
8tt 



pM , pM , pM , ^_ p , p +p 



which can be rewritten as 



B 2 

A\i e n e + AfipUp + Afi n n n + — = 

87T 



(47) 



(48) 



when the magnetic pressure is a small fraction of the total pressure. (As before, the 
superscript M denotes the interior of the flux rope, and A/x e = fj,^ 1 — /i e , etc.) The condition 
of charge neutrality, — np 11 = n e — n p = 0, leads to the relation 

UeAfle _ n p Afi p 
^ — i 

H e Hp 

and the condition of neutral buoyancy, n^f + n 1 ^ = n n + n p , leads to 
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n n Aii n UpA/jp 

ftn ftp 



= 0. 



(50) 



One solves for A/z by combining these three equations, 



A(JL 



2 fin 



(1 - Y e )fi ( 



- 1 



B 2 



8irn e 



(51) 



This expression can be simplified further after observing that the shift in the charged 
particle density is related to A/i via 



A/i = — 

3n e 



2ftn 



- 1 



so that 



M 

lip ~ Up 



1 - Ye) ft 

AY e 3 B 2 



rip -rip), 



4 8nP P 



(52) 



(53) 



This shift in the electron fraction is similar, in order of magnitude, to that induced by the 
Lorentz force acting on a homogeneous n — p — e plasma without gravity (Goldreich and 
Reisenegger 1992). 

The vertical speed Vr of the flux tube is limited by the rate at which the charged 
particle density is driven to its /^-equilibrium value by absorption of v e and v e on free 
nucleons, 

&Y e T r dY e T(v e + n -> e~ + p) - T(u e + p -> e+ + n) 



dt 



R 



OR 



p/m n 



(54) 
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Here T is the rate per unit volume, so that 

T(v e + n—>e~+p)- T{v e + p —> e + + n) = {a Uen )n v n^ c - (a DeP )n Pe n^ c. (55) 



In equilibrium, 



and so 



\&v e n)Tlu e Tip _ Y e 



\Vv eP )n Ve n n 1-Y e 



(56) 



T(u e + n —> e + p) - T(v e + p -> e + + n) = (cr Uen )n Ue (n™ - n n )c - (o DeP )n De (nf - n p )c 

= ^ 1 KnKc« -Tin). (57) 

Since AY e /Y e is assumed to be small, one can express Vr in terms of the gradient of the 
unperturbed electron fraction, 

T(v e + n -» e~ + p) - r(P e + p -» e+ + n) 
R ~ (dY e /dR)p/m n ■ 1 j 



Combining equations (p3f), (|57| ) with (p8|), one obtains 



Vr 3 11 ' ! 



c AY P 



d\nY e 



d\nR 



(<T Ven )n, e R) -f-. (59) 



Note that (a Ven )n Ue R = {n Ve /n n )r Ve = 7x 10~ 4 T Ue (T Me v/20) 3 (n n /n sat ) -1 . In this expression, 
r„ e = n n (a Uen )R ~ 10 4 (k B T/20 MeV) 3 is the absorption depth at 15 MeV £ A; B T ^ 40 
MeV (cf. Fig. 7 of Reddy et al. 1998). Substituting these expressions into (0), one obtains 

t" 4x10 Mtti5*gJ hr) (oos) ■ (60) 

At the Kelvin time of ~ 3 s, the density is n n ~ n sat and the temperature is k B T ~ 20 MeV 
at an enclosed mass of ~ 1 M & (Fig. 9 of pons et al. 1999|) . A flux rope is able to overcome 
the stable stratification from this depth only if V B /c ^ ^R QS /ct ~ 3 x 10~ 6 , where _R ns ~ 10 
km is the stellar radius. The lower bound on the flux density which can rise buoyantly is 
then 

5^3 xio 13 (%y 3 (^y 1/2 g. (6i) 



20 J V3 s, 

Notice the strong dependence on T. The direct URCA reactions freeze out at temperatures 
below k B T ~ 10 MeV, and so transport of the magnetic field to the stellar surface becomes 
ineffective at t ^ 10 s. 

The net effect of this transport process is to exchange a region of high flux density with 
a region of lower flux density that lies above it. Over such short timescales, the neutron 
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matter moves with the magnetic field. A field that uniformly threads the neutron core 
does not induce buoyancy forces and can exceed (|61|). Perhaps coincidentally, this bound is 
comparable to the dipole magnetic field inferred for the Soft Gamma Repeater sources in 
the magnetar model (Thompson and Duncan 1996; Kouveliotou et al. 1999| ). 



4. Magnetostatic Equilibria 

The convective core of a neutron star is surrounded by a convectively stable layer that 
extends from large to small neutrino optical depth, and straddles the neutrinospheres (Fig. 
1). A rope of magnetic flux that is forced from below by convective overshoot into this layer 
can continue to rise buoyantly above the neutrinosphere, through the effects of neutrino 
heating as discussed in §^|. We now consider the equilibrium shape of such a rope when 
its two endpoints are pinned from below (Fig. 3). The picture of the magnetic field that 
we adopt is one of spatially intermittent fibrils that are confined near the boundaries of 
convective cells below the //-sphere. Such a distribution is observed at the Solar photosphere 
( BtenHo 1989| ) and is inferred in the upper convection zone from helioseismology ( poldreicFi 



|et al. 1991| ) . It is also observed in simulations of a turbulent, conducting fluid at high 
magnetic Reynolds number ( [1'helen fc Cattaneo 2 



Aside from its dramatically higher temperature and stronger gravity, the exterior of 
newly formed neutron star is distinguished from the non-thermal atmosphere of the Sun in 
an important respect. Before the shock succeeds, the pressure remains high enough that 
magnetic flux originating inside the star will be confined to narrow fibrils well outside the 
//-sphere. To see this, we compare the pressure in the r.m.s. magnetic field flip]) inside the 
//-sphere, with the pressure of free nucleons near the //-sphere, 

(B 2 ) R j8rc fe B \ ( p \ X (T^y 1 ( R v \~ 3 

{P) = nk B T = -° 05 lolj (lOi'gcm-3 J UkJ " (62) 

This low volume average of /3 originates in the relatively low convective Mach number of 
the supernova core (as compared with the upper convection zone of the Sun), combined 
with the assumption of a common dynamo efficiency B 2 /4npV 2 on . 

The field that pierces the //-sphere is not force free. We now calculate the equilibrium 
shape of a flux rope that is underdense compared with its surroundings and confined by the 
external pressure, so that B 2 /8tt ~ P. The diameter of such a confined bundle of flux (of 
cross-sectional area ira 2 ) increases slowly with distance R from the center of the neutron 
core, 

4 oc R^~ l ~R°- rV\ (63) 
R 
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since the pressure stratification outside the //-sphere is P(R) oc R 13 {(3 ~ 4 — 5). 

As a result, the pressure scale height is a much larger fraction of the radius, £p/R ~ ~, 
near the v-sphere of a nascent neutron star, than it is near the Sun's photosphere, where 
£p/Rq ~ 10 -4 . As a result, a magnetic flux rope is subject to a strong buoyancy force on 
both sides of the ^-sphere. To calculate the equilibrium shape of the rope, we first consider 
the case of a plane-parallel atmosphere. A flux rope sitting in the x — z plane (Fig. 3), with 
tangent vector 

1 = cos 6*x + sin 9z, (64) 
is subject to a buoyancy force density 

(Pm ~ p)g = PpP\q\ r, (65) 

where 

d2 

/3p = — — — (66) 



and pm is the density inside the flux rope. In equilibrium, the normal component of (|65|), 

(Pa/ - P) [g - (g ■ 1)1] , (67) 
is balanced against the (normal) tension force density 

When neutrino heating equilibrates the temperatures inside and outside the flux rope 
(§|3.2|), the trajectory of the rope obeys the simple equation, 

d l = -JE- cose= -^l(P). (69) 

dl 2i m Zip \r m J 

Here P is the total pressure, and P m = pT/m n is the pressure in free nucleons (assuming 
that T > 1 MeV and alpha particles are absent). A related expression has been written 
down by Parker (1979) for a plane-parallel atmosphere with negligible pressure in relativistic 
particles. The main consequence of the pairs and photons is to reduce the scale height from 
t P = P/pg to 

Ip = y g = &) lp - (70) 
The separation between the two footpoints is 



Ax = J cosM = 2 J ^ £ P [l(9)]de, (71) 
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and the maximum height attained by the arcade is 

2 m ax = 2 tan elp[l(6)}de, (72) 
Jo 

where 8q is the inclination of the flux tube from the horizontal at z — 0. In an atmosphere 
with constant pressure scale height £p, these expressions simplify to 

Ax = 2l P 6 < 2ttI p , (73) 

and 

z max = 2l P \n( -). (74) 



\ COS 9o / 

The flux arcade expands to infinity as the footpoint separation approaches the critical value 

Ax c = 2nI P , (75) 

at which the flux segments become vertical, 9 — | . A similar effect will occur in any stably 
stratified, plane parallel atmosphere. 

Notice that radiation pressure can have a significant effect on the equilibrium shape of 
a magnetic arcade, since Ip oc pT/P. As a result, the critical footpoint separation decreases 
as photons and pairs contribute more to the pressure. 

These results are easily generalized to an atmosphere with spherical symmetry. The 
shape of a flux arcade that sits in the equatorial plane (with azimuthal angle 0) is defined 
by the equation 

di 2£ P 

which simplifies to 

# \2£ P J 

In a spherical atmosphere with I P /R a constant, the flux arcade attains a maximum radius 

I \ (R/2?p-l) _1 



-Rmax — Rv 7T , (78) 

VCOS^o/ 

and opens to infinity when the angular separation between the footpoints approaches the 
critical value 

A0 C = 2 S = —i . (79) 

Rj2i P -\ RjUp-X 
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For example, if the pressure is dominated by relativistic particles then lp ~ R/4 and lp is 
determined by the post-shock entropy. Substituting l P ~ R/5 outside the i/-sphere (around 
100 msec after bounce), one deduces a critical separation 




Magnetic arcades that extend beyond the i/-sphere are not fixed in position: their 
footpoints can be expected to wander in response to the convective motions below. Since 
these footpoint motions occur on the short timecale ~ r con , it is important to check whether 
an individual arcade is able to reach a new equilibrium configuration over the timescale 
~ Tcon- Otherwise, the magnetostatic approximation is not justified. The motion of the 
upper part of an arcade is driven by the tension force /tension = (na 2 )(B 2 /SnR c ) (where 
R c is the curvature radius and 7ia 2 the cross-section of the flux rope) and is resisted by 
the hydrodynamical drag force /drag = 2CdapV 2 . Given a footpoint separation ~ 2R C , the 
velocity of the rope (unimpeded by drag) is V ~ 2R c /r con , and the ratio of drag to tension 
forces is 

/drag = QQ J_ ( R \ 2 [TueVy 1 f2R c \ 2 ( T con \ ~ 2 / a 

/tension ' Pp V50 km) V 4 J \ i P J \3 ms) \0.5R C 

One sees that smaller flux arcades (2R C ~ lp) that sit near the //-sphere [R v ~ 30 — 50 km) 
will be able to follow the footpoint motions, whereas larger arcades that extend to a radius 
R^> R v will lag behind, and probably become very tangled. 

Self-reconnection is also possible for smaller flux arcades that expand into the 
convective layer that sits below the shock; or for those which rotate rigidly with the central 
neutron star and experience a ram pressure force density 

p(R)(n coie R) 2 (82) 

that exceeds the buoyancy force density ( |65|) . This second condition is satisfied beyond a 
radius 

*- = 320 (TS^) 2/3 fe) 1/3 (^) V3 < 83 » 

where P cor c = 27r/f2 core is the rotation period of the neutron core. However, as discussed in 
Chevalier (1989) and Thompson (2000), the neutron star may acquire its rotational angular 
momentum only during the last stages of accretion from the pre-supernova core. 
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5. Implications for Pulsars 

This paper has investigated the response of a magnetic field to the violent fluid motions 
and intense neutrino flux within a supernova core. We can summarize our conclusions as 
follows: 

1. Two principal sources of free energy are available for amplifying a magnetic field in 
a non-rotating supernova core: the violent convection occuring within the neutrinosphere 
(TD93); and within the gain region in between the neutrinosphere and shock ( [100| ) . The 
equilibrium magnetic field in these two convective zones is respectively ~ 10 15 and ~ 10 14 
G (eqs. |TJj and [fTTH). The convection occuring within the neutrinosphere is a robust 



consequence of the diffusive transport of heat by neutrinos (Appendix A). Even during 
relatively late fall-back (~ 0.1 M over ~ 10 3 s), the neutrino flux is high enough to induce 
vigorous convection within the material below the accretion shock. 

2. The supernova core also contains regions which are stable to Ledoux convection, 
including the material which straddles the neutrinosphere(s). The transport of magnetic 
fields across these regions is greatly facilitated by the absorption of electron-type neutrinos 
(and anti- neutrinos) on free neutrons (and protons). The field induces a chemical potential 
imbalance between electrons, neutrons, and protons (cf. Goldreich & Reisenegger 1992), 
but the charged-current reactions force this imbalance toward zero (§§). During the main 
Kelvin- Helmholtz cooling phase of a nascent neutron star, this process effectively transports 
material containing fields stronger than ~ 10 13 G out of the star. 

3. During periods of hyper-Eddington accretion, the equililibrium shape of a magnetic 
flux rope which protrudes across the neutrinosphere is strongly modified by the hydrostatic 
pressure of the settling material. In the case of a flux rope anchored at both ends, the 
resulting buoyancy force leads to a critical separation between the ends of the flux rope, 
above which the rope expands to very large radius (§§). This mechanism provides a means of 
feeding large amounts of magnetic flux from a dynamo operating within the neutrinosphere, 
to the neutrino-heated 'gain region' outside the neutron core. The interaction of this field 
with the accretion flow has been investigated elsewhere ( |TO0| ). 

Claims that late fallback will tend to bury a magnetic field (e.g. Geppert, Page, & 
Zannias 1999) should, as a result, be treated with caution. The full implications of these 
physical processes for pulsar magnetism (and perhaps also the supernova mechanism) can 
only be gleaned by incorporating them into full numerical simulations of the supernova 
collapse. The behavior of a centrifugally-supported accretion flow is, in particular, beyond 
the scope of this paper. Nonetheless, several qualitatitive conclusions will be summarized 
here, in the following sections. 
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5.1. Principal Phase of Dynamo Amplification 

Both the gain region below the accretion shock and the interior of the neutron star 
are strongly convective. The ratio of convective kinetic energy to gravitational binding 
energy is a hundred times larger in the proto-neutron star phase, than during any previous 
convective phase driven by nuclear burning (TD93). It is even larger in the convective gain 
region below an accretion shock flT00| ) . Since the ratio of magnetic energy to gravitational 



binding energy is approximately conserved during gravitational collapse, these violent fluid 
motions during the supernova event are expected to be the dominant source of free energy 
for pulsar magnetism. The relative importance of the two convective regions for the surface 
fields of radio pulsars, depends in part on the amount of fallback following the success of 
the supernova shock. 

Even during late infall (10 3 — 10 4 s after collapse) the settling flow is heated 



sufficiently to become violently convective. As reviewed in § [2.2| , the settling speed scales 
as V(R) oc M 3 / 5 when cooling is dominated by e ± captures on free protons and neutrons, 
whereas the convective speed decreases more slowly, V con (R) oc M 1 / 3 (at fixed radius R). 
The equipartition magnetic field within the flow scales as (B 2 ) 1 ^ 2 oc M 8//15 (eq. |p0[| ). 



5.2. Prompt Convective Dynamo and the Surface Magnetic Fields of Young 

Neutron Stars 

The bulk of the crust of a neutron star is derived from material that is either processed 
through the convective shell below the accretion shock; or rises buoyantly from the deep 
interior of the neutron core through the convectively stable layer that straddles the i/-sphere. 
The equilibrium flux density (|T7|) in the outer convective shell, scaled to the density of the 
neutron star crust, is 

B crust ~ I (5 2 ) 1/2 , (84) 

or equivalently, 

9 ln 17f ( M V 1/6 / R sh V ^ 2 fC^' 2 ( M COT e p 

(85) 
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Here p 2 is the post shock density (eq. flIE|). This field can be expected to be coherent on 
the small scale 

f ( V /3 t cr \ n , ( R ^ V V2 ( & ^ ( n \ 1/3 ( M ™ v 

' fl ~UJ ^ (JW = -HiootoJ \Tm^) (ioJ (UmoJ km - 

(86) 

This leads to a net dipole field 

B«p* ~ <5 2 >V 2 (-^j , (87) 

namely 

(88) 

A further multiplicative factor of (Pb/Rns) 1 ^ 2 ~ (P2/Psat)^ 6 must be included if one also 
averages the dipole field through many radial shells of thickness £b- This estimate applies 
to the field generated by convective motions in the outer part of the gain region. 

Material that rises buoyantly past the v-sphere must also be strongly magnetized. The 
strong neutrino flux aids buoyancy forces, but the magnetized fluid will reach the surface of 
the star during the prompt Kevlin-Helmholtz phase, only if B > 10 13 G (eq. pl|). Weaker 
magnetic fields that are more characteristic of the dipole fields of ordinary radio pulsars 
are not able to reach the neutrinosphere, and according to our current understanding of 
transport through the crust will remain buried for ~ 10 8 yr or longer (e.g. Goldreich 
& Reisenegger 1992). This suggests that a prompt convective dynamo can explain the 
observed magnetic dipole fields of radio pulsars only if the dynamo operates stochastically - 
generating a true dipole component through the incoherent superposition of stronger, small 
scale fields (TD93). 



5.3. Effects of Late Fallback 

Even after the shock is pushed to large radius, continuing fallback (driven by the 
formation of a reverse shock) can keep the hydrostatic pressure at the surface of the nascent 
neutron star well above the surface value of B 2 /8n inferred for the dipole fields of young 
radio pulsars. For example, the hydrostatic pressure of a spherical shell of mass AM 
corresponds to a field strength 

/ 2GMnsAAA ' /2 , lnl7 / AlV" I M m V /2 ( Rm \~ 2 „ ,„„ 

B= { w s ) =3xl ° [oiwj [tam;J iw^> g - < 89) 
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We focus here on the consequences of spherical accretion. (The behavior of a centrifugally 
supported flow may be substantially different, because it may not be able to maintain the 
large central cusp in the temperature and density needed to effect rapid neutrino cooling.) 

Consider first the equipartition magnetic field which is maintained by the convective 
motions within the material settling onto the neutron star, driven by neutrino heating. At 
low accretion rates, the accretion shock lies at a large radius, and so it is appropriate to 
focus on a fixed radius R 100 km, where (B 2 ) 1 ^ 2 oc M 8//15 (eq. |2T)fl). Allowing for an 
additional factor of ~ (psat/p) 1 ^ 6 - relating the net dipole field to the equipartition field in 
the settling flow (§ |5.2| ) - one obtains the relation 



B dipole oc M 7/15 (90) 

between the dipole field generated by convective fallback, and the rate of mass accretion. 
The estimate fl88|) then reduces to -Bdipole ~ 10 13 G at an accretion rate M = 1O _4 M s _1 . 

Next, we consider how the continuing flux of electron- type neutrinos will promote the 
buoyancy of magnetic fields which are anchored below the final layer of accreted material. 
We start with eq. (J33|) for the radial drift speed of a (horizontal) magnetic flux rope through 
an atmosphere of non-degenerate nucleons, given a heating rate ( P7D by charged- current 
absorption of v e and u e . At an intermediate accretion rate 10 _1 ^ M ^ 10 -4 , the 
gravitational binding energy is carried off through e ± capture on free nucleons in the 



settling flow (|T00|) , as compared with e annihilation at lower accretion rates (Chevalier 



1989). The associated luminosity is directly related to M by eq. (0). Making use of the 
scalings ( |18"D of density and temperature with M (at fixed radius within the settling flow), 
one finds that the nucleons become less degenerate as M decreases. Magnetic flux poking 
out of the nascent neutron star will, as the result of buoyancy forces, carry with it a slightly 
lower mass density and maintain a similar degree of degeneracy to the surrounding flow. 

It is now straightforward to derive the critical (3p = B 2 /8imT above which a parcel of 
magnetized fluid will rise buoyantly over the infall time Ti n f a u ~ 10 3 s. One obtains 



a,* i. 6 2=! =ixio- ( 4 m y i/iQ f^r i 

H P< Tinfoil V 10 M © 8 / U0 3 S y UOW V L4M ©/ 

from equations (|27|) and (|43|) . The tiny numerical factor shows that the surface magnetic 
field almost certainly cannot be shielded diamagnetically by the accreted material 
(AM ~ Mr in f all ), because in that circumstance (3p ^ 1. Burial is possible only if the 
pre-existing magnetic field undergoes partial turbulent mixing with a fraction of the 
accreted material, so as to force [3p below the value (EH1). 
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5.4. Origin of the Low-order Magnetic Multipoles of Neutron Stars 

At least two distinct types of hydromagnetic dynamos can give rise to the dipole 
magnetic fields of neutron stars. The dipole field can either be amplified by fluid stresses 
that are coherent over a substantial fraction of the stellar radius (e.g. the classical a — Q 
dynamo driven by a combination of rapid differential rotation and convection); or it can 
be built up by much smaller units that themselves are more strongly magnetized. In the 
convective interior of a neutron star, the first type of dynamo requires very rapid rotation, 
P ^ ^p/V con ~ 3 ms, and has been conjectured to result in magnetars with dipole fields in 
excess of 10 14 G. The second type of stochastic dynamo has been conjectured to generate the 
small scale, intranetwork magnetic field of the Sun (TD93; |Durney, de Young, fc Roxburgh 



1993| ), but is ineffective at generating low order magnetic multipoles in a main sequence 
star. For example, the pressure scale height at the Solar photosphere is only ~ 10~ 4 of 
the Solar radius, which implies that the dipole flux density resulting from the incoherent 
superposition of the small scale dipoles is suppressed by a factor ~ (2nR 2 & / l P )~ l l 2 ~ 10~ 4 . 
By contrast, lp ~ i? ns /30 at the neutrinosphere of a newly formed neutron star, during 
the last stages of Kelvin-Helmholtz cooling. The corresponding suppression factor is 
~ 10~ 2 . This can comfortably accommodate the dipole fields of young radio pulsars, which 
characteristically lie in the range ~ 3 x 10 11 — 10 13 G. 

A similar stochastic dynamo may operate in the convective accretion flow onto a 
nascent neutron star. During •prompt fallback, within ~ 1 s after collapse, this mechanism 
could generate a dipole field as large as ~ 10 14 — 10 15 Gauss (eq. H88|). However, the 



equipartition magnetic field in the convective settling flow becomes weaker with increasing 
fallback time (and decreasing M; eq. This means that variability in the amount of 

fallback is an alternative explanation for the presence of a range of dipole fields in young 
neutron stars - with magnetars representing those objects with the smallest amount of late 
fallback. It is interesting to note that the current (uncertain) estimates of the mass accreted 
during late fallback (AM ~ 0.03 - 0.1 M ; e.g. Fryer, Colgate, & Pinto 1999) typically 
exceed the mass of the rigid crust in model neutron stars (~ 0.02 M ). As a result, the field 
generated during fallback will tend to be anchored in the neutron star core. 



5.5. Observational Tests of the Dynamo Mechanism 

Given the extreme difficulty of inferring the initial spins of radio pulsars, models of 
the dynamo origin of pulsar magnetic fields are most tightly constrained by searches for 
correlations between B and O. The stochastic dynamo does not depend in any way on 
rotation, and so predicts that these two vectors are randomly oriented with respect to 
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each other. By contrast, the magnetic field formed by a large-scale helical dynamo in 
a rapidly rotating neutron star should be strongly correlated with the rotation axis. In 
addition, a neutron star containing a toroidal magnetic field stronger than ~ 10 14 G (and 
an external magnetic dipole that is approximately aligned with the symmetry axis of the 
internal field) has a tendency to reach an alignment between Q and the external magnetic 
moment (cf. Goldreich 1971). This effect is more likely to induce alignment in the Soft 
Gamma Repeaters and Anomalous X-ray Pulsars, than in young radio pulsars (Duncan, Li 
& Thompson 1994). 

5.6. Relaxation of Small-Scale Magnetic Fields 

A small scale magnetic field of strength 10 14 — 10 15 G in radio pulsars would be 
manifested through sudden fractures of the neutron star crust. This phenomenon could 
possibly play a role in triggering glitches (TD93) by intermittently heating the crust, 
or inducing rapid translations of the Coulomb lattice with with respect to the neutron 
superfluid. Nonetheless, the absence of X-ray bursts associated with radio pulsars, combined 
with the seismic quiescence of old radio pulsars, suggests that if such a field is present 
it must be buried fairly deeply. Either a radio pulsar must have accreted upwards of 
~ 10 -2 M of weakly magnetized material (B ~ 10 12 — 10 13 G) following the brief dynamo 
epoch (eq. [pOfl), or the high order magnetic multipoles produced by convection were 
effectively smoothed out. 

The stability of higher magnetic multipoles depends on the ability of the gravitationally 
stratified fluid to undergo fully three-dimensional motions (TD93). Stable stratification, 
which forces the fluid to move along the equipotential surfaces of the star, may allow 
more complicated magnetic topologies. However, we have shown that the absorption and 
emission of electron-type neutrinos facilitates the buoyant motion of magnetic fields as weak 
as ~ 10 13 — 10 14 G (during the prompt Kelvin- Helmholtz phase of a proto-neutron star), 
and so will facilitate the reconnection and unwinding of stronger fields. 

ACKNOWLEDGEMENTS 

We gratefully acknowledge conversations with Phil Arras, Leon Mestel, Ewald Mueller, 
and Mai Ruderman. Financial support has been provided to both authors by the NSERC 
of Canada; and to C.T. by the Sloan Foundation. 



-32 - 



A. Convection in Young Neutron Stars 

In this appendix we calculate the temperature gradient required to trigger convection 
in the presence of the stabilizing composition gradient found in a cooling neutron star. 
This gradient is derived below and given in eqn. (H) in the main text. We compare the 
result to that of Lattimer & Mazurek (1981), who find that lepton- fraction gradients 
tend to be stablizing at low lepton fraction, and Reisenegger & Goldreich (1992), who 
find that the Y" e -gradient in a cold neutron star at /9-equilibrium is absolutely stabilizing. 
The critical value of Yi, below which the lepton gradient is stabilizing, increases at high 
temperatures and low densities, a result embodied in eqn. ( |A33[ ) below. Although more 
detailed numerical cooling models are now available (e.g. Pons et al. 1999 ), the following 
arguments are useful for the physical insight they provide. 

We are interested in the response of a neutron star to displacements of fluid elements. 
The temperature and density of the stellar material vary with radius; moving material 
radially will result in buoyancy forces which can either oppose or enhance the force initiating 
the displacement, depending on the thermodynamic state of the star. We assume that the 
star is in hydrostatic equilibrium, with equilibrium density and pressure gradients dpo/dR 
and dP /dR, related by 

-Po9 (Al) 



dR 

and an equation of state 

P = P(j>,T,Y e ,Y Ve ). (A2) 

The density and temperature are assumed to be high enough that the stellar material 
is composed entirely of neutrons, protons, electrons, and neutrinos.^ The relative 
proportions of these species are related to the total nucleon density^] n& = p/m n by the 
parameters Y e = (n e - — n e +)/nb, Y p = n p /nb = Y e (from overall electric charge neutrality), 
Y Ue = {n Ve — n Pe )/n b , and Y n = n n /rib = 1 — Y e . The total lepton fraction is 

Y l = Y e + Y Ue . (A3) 

We model each species of particle as an almost degenerate, ideal fermi gas with chemical 
potential pi. The background configuration is assumed to be in chemical equilibrium: 

Pn + Pu e = Pp + p e - (A4) 



4 But not so high that p, and r leptons, pions, or strange particles are generated in significant abundances. 
5 We ignore the electron mass as well as the small mass difference m n — m p . 
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In a hot neutron star (T ^ 10 MeV), which is optically thick to all the neutrino species, 
a displaced fluid element can be assumed to remain in chemical equilibrium. At very 
high temperatures, the lepton fraction can also be assume to be frozen into a buoyant 
plume of material, Y\ ~ constant. The convective efficiency is significantly reduced once 
the optical depth r Ve to electron neutrinos across a pressure scale height drops to a value 
~ c/Kon ~ 300, where V con ~ 10 8 cm s _1 is the speed of the convective overturns (cf. 



Mezzacappa et al. 1998|) . The corresponding mean free path to v e + n — > p + e~ is A ~ 10 3 
cm, which is reached at a temperature ksT ~ 15 MeV at nuclear saturation density 
(rib = n sat = 1-6 x 10 38 cm -3 ; cf. Fig. 7 of Reddy et al. 1998). In fact, as we now show, this 
temperature is comparable to the value below which a purely radiative flux of neutrinos can 
be maintained without inducing convective instability. 



A.l. Convective Instability 



We first review the derivation of the buoyancy force and the Brunt- Vaisala frequency 
N associated with the displacement of a fluid element from its equilibrium position. After 
a displacement through a distance £ in the radial direction, the density of a fluid element 
becomes 

dp "N 



Pd (R + = p(R) + 



dR, 



(A5) 



S,Yi 



The subscripts indicate that the entropy and composition are held fixed in calculating the 
derivative. 

The buoyancy force arises because this density differs from the background density at 
r + £: 

p (R + C)=p(R) + ^. 



The (Eulerian) density difference is 



8p = p d (R + ^)-p (R + C)^ (J£l 



dR 



dp 



dR 



resulting in a buoyancy force 

F B = -g5p = -g 



dp " 
dR 



S,Y, 



S,Y, 



dpp 
dR 



dR 



(A6) 



(A7) 



(A* 



Suppose we move a bit of fluid upward, so that £ > 0. If the density of the displaced fluid 
element is less than that of the surrounding material (5p < 0) then the term inside the 
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square brackets in eqn. ( |A8|) is negative. The buoyancy force is upward, in the direction of 
the original displacement, and the star is unstable to convection. 

Young neutron stars have both entropy and composition gradients, so it is convenient 
to express the buoyancy force and N in terms of those gradients. Noting that 
p(R) = p [P(R), S(R), Yi(R)] we find 



' dp N 
dR, 



S,Yi 



and 



dp " 
dP 



S,Y 



dPo 
dR 



dpo 
dR 



[dp / 
'dp) 

dS 



S,Y, 



dPo 
dR 



P,Y 



dS 
~dR 



+ 



' dp " 
W 



p,s 



dYifl 
dR 



(A9) 



(A10) 



The condition for convective instability becomes 



' dp s 
dS 



P,Yi 



dS 
~dR 



' dp \ dY lfi 
dYj dR 



> 



Since y§gj py < 0, the condition for triggering convection is 



dSo (dp/dY^s dY lfi = dSo_ ( dS_\ 
dR (dp/dS) P , Yl dR dR [dYj 



dYifl 
dR 



< 0, 



(All) 



(A12) 



which is eqn. ([!]) in the main text. 



A. 2. Influence of a Compositional Gradient 

The sign of the thermodynamic derivative (dS/dYi)p tP determines whether the lepton 
gradient is stabilizing or de-stabilizing. Note that dYi/dR < is quickly established by 
passive neutrino diffusion in a proto-neutron star. Evaluation of the derivative is greatly 
simplified when Y v <^ Y e , so that one can approximate d/dYi ~ d/dY e . 

We idealize the neutrons and protons as non-relativistic Fermi gases. It will prove 
convenient to introduce a characteristic Fermi momentum 

1 /3 

p F = h(3n 2 n b ) 1/3 = 330 ( — ) MeV/c. (A13) 

V^sat/ 

Then the Fermi energies of the neutrons and protons can be written as 




(A14) 
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Here % is Planck's constant. The neutrons are in fact just becoming degenerate at k B T « 40 
MeV, which is the peak temperature generated immediately after the collapse ( [Pons et~al 



1999[); it is the associated rise in pressure that halts the collapse. The protons are less dense, 



and so they are only marginally degenerate. Their contribution to the entropy is smaller 
but significant. The analysis is simplified by treating the protons as an almost degenerate 
gas, but it should be kept in mind that this approximation is only marginally valid. 

The electrons and electron neutrinos are, by contrast, highly relativistic (cf. equation 
A13|| ). Their Fermi energies are 



e e « hc(3n 2 n e ) 1/3 = F e 1/3 p F c. 



and 



Notice that 



ftc(67rV) 1/3 = (ZnjVs 



p F c. 



~y 



(A15) 
(A16) 
(A17) 



is not necessarily small when Y Ve <C Y e . We will not neglect terms of order e Ue /e e even while 
dropping terms of order Y Ue /Y e . 



At finite temperature, the chemical potentials are 



and 



7T 

12 



3 



k B T 



and the entropies (per baryon) are 



S 



p,n 



_ J T 



p,n 



and 



e , v 



k B 

The pressures can be approximated as 



7r 2 y fi . 



-p,n 



k B T 



and 



p,n ^^p,n^p,n r^bSp^-nT 



1 + 1 S^T. 

4 o 



(A18) 
(A19) 

(A20) 
(A21) 

(A22) 
(A23) 
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The key thermodynamic derivative can be expanded as 



'dS\ 



p,p 



'dS 

W P 



l p,p \ dY ^ 



+ 



or 



P,Ye 



dT 
dY 



p,p 



To evaluate this expression, we need 



'dS_ s 
dT 



P,Ye 



S_ 

T' 



and 



'dT 



(dP/dY e ) p , T 

] Pp - (dP/dT) pXe - 



(A24) 



(A25) 



(A26) 



The derivative (dP/dY e ) TjP is dominated by the (zero-temperature) degeneracy pressure 

2 2 11 

= -n n t n + ~n p e p + -n e e e + -n Ue e Ue , (A27) 

so that 



1 , ,1 dY u 

-n b {2e p - 2e n + e e ) + ~n b e Ue I 



(A28) 



Neglecting the last term on the RHS for the moment, and invoking the /3-equilibrium 
condition e p + e e = e n + e u , one finds 



(dP \ n b 

However, the thermal contribution to the pressure cannot be neglected, 



(A29) 



'dP_ 

dY P 



P,T 



'dPo 
dY„ 



+ 



n b T d 
~^~dY 



S p + S n + — (5 e + S v 



(A30) 



as it rescales the first term on the RHS of eq. ( |A24j ) . Combining 



'dP s 
dT 



P,Ye 



n b [ S p + S n + -S e + -S u 



(A31) 



with the above equations, one finds 



'as 

dY 



3 
+ 4 



'OSe 

dY, 



+ 



P,T 



P,P 

( dS Ve 
\ dY 



'dS n 



dY P 



+ 



P,T 



'dSp 
dY P 



PP. 



S p + S n 



l S p + S n + (S e + S 1/e )/2 



+ 



P,T 



S p + S n + (S e + S Ue )/3 
S p + S n + (S e + S Ue )/2 



S p + S n + (S e + S„ e )/2 



2e Ve - e e 
2T 
(A32) 
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Neglecting the entropy of the electrons and neutrinos compared with the nucleons, this 
expression simplifies to 



dS 



OS 



2e»„ - e e 



BY / V W 7 

, UIe J p, P " \ w±e / p,T 



2 \dYj „. 27 

The variation of the entropy with lepton number is, in the same approximation, 



(A33) 



OS \ 1 fS, 



dY e ) pT 3\Y e l-V; 




1 - 



Y, 



l-Y P 



2/3' 



(A34) 



A fuller evaluation of these derivatives, including the contributions from the 
electrons and electron neutrinos, requires knowing the quantity (dY Ue / dY e ) pT , 
which can be determined directly from the condition ( |A4j ) of chemical equilibrium: 
(dY Ve /dY e ) p> T — (Y„ e /Y e ) 2 ^ 3 when Y Ve <C Y e . For example, substitution of this expression on 
the RHS of eq. ( |A28|) shows that the derivative of the electron neutrino pressure is of order 
Y Ue /Y e , and can indeed be neglected. 



The compositional gradient (dYi/dR < 0) tends to suppress convective instability if 
(dS/dY[) PtP > 0, and is otherwise de-stabilizing. Both terms in eq. ( |A33| ) 
2e v < e e , that is, if Y Ue < jqY £ . However, immediately after the collapse e 



are positive if 
, ~ e P , and so 



the composition gradient is potentially de-stabilizing. The sign of (8S/dYi)p :P reverses at a 
critical temperature 



which we evaluate in the main text. 







1/2 


/ v \ 2//3 " 


-1/2 


7T 


m n c 2 
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This expression qualitatively reproduces the domain of convective instability in the 
S — Yi plane, as plotted in Figure 1 of Lattimer & Mazuerk (1981). Those authors find that 
low Y\ and high T tend to suppress compositionally-driven convection. They also find that 
convective instability remains possible at very low entropies and lepton fractions, where 
heavy nuclei can form. However, that instability is restricted to densities below nuclear 
saturation. 



A. 3. Triggering Convection 



When the compositional gradient is destabilizing, convection will ensue. However, when 
it is stabilizing, convection may still occur if the entropy gradient is negative. Convection 
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occurs if 



dS 
~d~R 



< 



'd£T 



dYifi 
dR 



(A36) 



The gradient of YJ is generally negative, and we have already evaluated the 
thermodynamic derivative (dS/dY{)p t p. It remains to re-express the entropy gradient in 
terms of a convective temperature gradient. Approximating 

•-' V • 

(A37) 



S*S n + S p = 2%££ [yV3 + ( i_y e) i/3 
Sn rib 1 



(recall pp = h (Sir 2 rib) 1 ^ 3 ), the left hand side of eqn. ( |A36| ) yields 
d In S dhiT 2 d In p 
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dlnR 



(A38) 



We evaluate the first factor on the left hand side of eqn. (|A36[ ) at low temperatures. This 
allows us to neglect the first term in eq. |[A33j| , yielding 



y p fas 



3 y4/3 



h cn b ( 2fi Ue - fji e 



(A39) 



S \dY e ) Pp yV3 + (1 _ y e) i/ 3 m n {k B Tf \ 2// e ) " 
Using eqs. ( |A38|) and ( |A39|) in (|A36|) leads to expression (§) for the convective temperature 



derivative (d\nT /d\n R) conv . In § [2.1| , we calculate the radiative temperature gradient 
associated with passive neutrino transport, and infer the critical radiative flux that drives 
convection. 



A. 4. The Brunt- Vaisala frequency 

The electron fraction decreases outward through the degenerate core of a cold neutron 
star. That this compositional gradient tends to stabilize against convection was noted by 
Lattimer & Mazurek (1981) and demonstrated in detail by Reisenegger & Goldreich (1992) 
(in the approximation where the neutrons and protons are normal Fermi fluids). Here we 
reproduce the result of Reisenegger & Goldreich (1992). 



The equation of motion for our fluid element is 

d\ „ 



P 



dt 2 



B- 



This leads to the definition of the Brunt- Vaisala frequency N, 



N 2 = 9 - 



dp s 
dr 
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(A40) 



(A41) 
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We use the latter form in order to check our expressions for the partial derivatives by 
comparing to the result of Reisenegger & Goldreich (1992), who evaluate the Brunt in a cold 
neutron star. The star is assumed contain no neutrinos (Y Ve = e v = 0), to be in chemical 
equilibrium (e n ~ e e ), and to have zero temperature and entropy. Under those conditions 
the Brunt reduces to 



T 2 9 ( dp\ ( dS \ dY e 



" = Z[i*)„\w.) p ^- (A42) 



From the definition Y e = n e /rib we have 

dY e 1 dn e Y e drib 



(A43) 



dr rib dr rib dr 
From the equilibrium condition we find 

an 2 n = n e , (A44) 

where a = (h/2mc) 3 37r 2 is a numerical constant. From charge neutrality n e = n p , 
conservation of baryon number gives n n + n e = rib. Taking the derivative of this expression 
with respect to radius, and solving for the electron density gradient, we find 

^ =2Ke ^ + ( n) 2 . (A45) 
dr dr 

Combining these results in 

**L = Xl d J*L. (A46) 
dr po dr 

Finally we evaluate the variation of density with entropy at fixed pressure and 
composition. We have 



'dS\ fdS\ fdS\ {BP I dp) 



T.Yi 



\dpJ PXl \dp) TYi \dT) pYi (dP/dT) pXi 
From equations (|A20[ - |A23|) we find 
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ar , Y r A (A50) 
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(A51) 



where we neglect the leptonic contribution to the entropy. Combining these with eqn. 
( TO we find 

fdS\ _ f2S\ f5S\ P 

The second term on the right-hand side of this equation is much larger than the first: 
P/n b ST « (4/5vr 2 )(e„/A; 6 T) 2 >> 1. To lowest order in fc 6 T/e 6 

\ds) pYe - 2\eJ S~ e n Y n {Ab2} 
Using equations (|A3|), ( |A46|) , and ( [A52D 

" po V e n Y n 2Tj po rfr 2y n #' 1 ' 

where H = —p/(dpo/dr). This is the expression found by Reisenegger k, Goldreich (1992). 
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Fig. 1.— The collapsed core of a massive star develops a violent convective instability 
in two distinct regions: within the central, neutronized core which is optically thick to 
neutrinos; and outside the neutrino photosphere, within a spherical shell where heating by 
the charged- current absorption of v e and v e on free n, p is faster than cooling by captures of 
electron-positron pairs. The region straddling the z/-sphere is, typically, stable to convection. 



-44- 




Fig. 2.— A horizontal bundle of magnetic flux has a lower temperature than its 
surrroundings, Tm < T, at constant density and pressure. Electron-type neutrinos, 
originating outside the bundle, raise its temperature and allow it to move upward through 
a stabilizing composition (or entropy) gradient. 



-45 - 




Fig. 3. — A narrow bundle of magnetic flux reaches an equilibrium configuration which 
is a competition between buoyancy and tension forces. When the separation between the 
magnetic footpoints exceeds a critical value, the equilibrium bundle reaches to infinite height 
(radius). 



